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Studying the interaction between microquasar jets and their environments 
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In high-mass microquasars (HMMQ), strong interactions between jets and stellar winds at binary system scales could 
occur. In order to explore this possibility, we have performed numerical 2-dimensional simulations of jets crossing the 
dense stellar material to study how the jet will be affected by these interactions. We find that the jet head generates 
strong shocks in the wind. These shocks reduce the jet advance speed, and compress and heat up jet and wind material. 
In addition, strong recollimation shocks can occur where pressure balance between the jet side and the surrounding 
medium is reached. All this, altogether with jet bending, could lead to the destruction of jets with power < 10^^ erg/s. 
The conditions around the outfiow shocks would be convenient for accelerating particles up to ~ TeV energies. These 
accelerated particles could emit via synchrotron and inverse Compton (IC) scattering if they were leptons, and via 
hadronic processes in case they were hadrons. 
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1. Introduction 

The interaction between a jet and its environment 
allows the study of both the properties of the jet 
and its ambient medium. ^'^ Moreover, this interac- 
tion is important for the evolution of jets.^ In the 
case of extragalactic jets, it is known to be different 
in case of FRI and FRII sources; the collision and 
interaction of jets with galactic clouds and collective 
stellar winds is thought to be partly responsible for 
the jet mass loading and deceleration in the former,^ 
whereas the latter can propagate and interact with 
the medium up to Mpc scales. In the case of micro- 
quasars, not much is known about the interaction 
of jets with their environments,^'^ and only several 
sources show evidence of such an interaction (e.g. 
SS 433;^ XTE J1550-564;8 Cygnus X-l^). 

In HMMQ, the primary star suffers severe mass 
loss in the form of a supersonic wind,^^ which can em- 
bed the jet rendering their mutual interaction likely. 
Here we report on a work studying this interaction 
and its effects on the jet evolution and radiation. 



2. Numerical Simulations 

We have performed three numerical simulations of 
jets with cylindrical geometry and two with slab (i.e. 
planar^) geometry propagating through the wind re- 
gion. We simulate the injection of a jet, with the 
properties given in Table 1, in the ambient medium, 
i.e. the stellar wind. The wind parameters corre- 
spond to those of a HMMQ similar to Cygnus X-1 
or LS 5039. The density is 2.8 x 10"^^g/cm-^, for 
a star with dM^/dt = IQ-^Mo/yr at 3 x 10^^ cm 
from the jet base. For simplicity, at this stage the 
wind porosity (e.g. Ref. 11) is assumed to be neg- 
ligeable. Typically, the velocity of an O star wind is 
~ 2 X 10^ cm/s. Despite this, the ambient is consid- 
ered to be at rest in the cylindrical simulations, this 
fact being, in principle, justified because the velocity 
of the jet is about 100 times larger that the wind one. 
It allows an axisymmetric treatment of the problem 
though the real wind is moving and comes from one 
side, being thus intrinsically asymmetric. We note 
that this assumption is not valid for weak jets on the 
basis of the results of the slab geometry simulations. 
The internal energy of this medium is obtained from 



^The slab jet has the same properties as in the cylindrical case per distance unit in the direction perpendicular to the simulated 
plane, and it is infinite in such a direction. 
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the temperature assumed for the wind (10^ K), with 
a resulting pressure of 1.5 x 10~^ erg/cm~^. The adi- 
abatic index of both the ambient and the jet (super- 
sonic and weakly relativistic), is 5/3. The parameters 
in the jets are summarized in Table 1. 

Table 1. Parameters of the simulations 



Parameter 



Weak j. Mildj. Powerful]. 



Jet power (erg s-^) 3.0 x 10^4 iQ^e 3.0 x lO^^ 

Jet pressure (erg cm-^) 9.1 68 6.2 x 10^ 

Jet density (g cm-^) 2.2 x IQ-^^ 5.9 x IQ-^^ 1.8 x IQ-^^ 

Jet speed (cm s-i) 1.3 x lO^o 2.2x10^0 2.2x10^0 



The numerical simulations were performed us- 
ing a two-dimensional finite- difference code based 
in a high resolution shock-capturing scheme which 
solves the equations of relativistic hydrodynamics. 
This code is an upgrade of that described in Ref. 1. 
The numerical grid is formed by 320 cells in the ra- 
dial direction and 2400 cells in the axial direction, 
with physical dimensions 20x300 jet radii. Extended 
grids were added in both directions in order to bring 
the grid boundaries far from the region of interest. 
The jet is injected at a distance of 6 x 10^^ cm from 
the compact object, and it is given an initial ra- 
dius of one tenth of this distance. During the simu- 
lation, the jet covers a significant fraction of the typ- 
ical size of a HMMQ binary system (here we adopt 
^orb ^ 3x lO^^cm). The planar simulations were per- 
formed with the same numerical resolution (which 
implies doubling the grid radially) and a shorter ax- 
ial grid by a factor 2/3. 

3. Results 

Figs. 1 and 2 show the mild and powerful jets when 
reaching the end of the grid, respectively. A large 
pressure jump (> 10^) is generated at the bow shock 
in both cases. A reverse shock in the jet itself is also 
formed. The mild jet generates a high pressure co- 
coon, which keeps it collimated, whereas the powerful 
jet develops a strong standing (recollimation) shock, 
with a pressure jump of the order of ~ 10^, due to 
overpressure with respect to its cocoon. Nevertheless, 
the mild jet also shows weaker internal shocks with 
pressure jumps of the order of ~ 100. The simulation 
of a weak jet (not shown here) gives results which are 
very similar to those of the mild jet. 



Dimensional estimates^^ show that, in the 
jet/cocoon phase, standing shocks must appear in 
the binary region (~ 10^^ cm from the jet base), at 
jet-temperatures > 10^^ K at injection in the nu- 
merical grid. The temperature of the mild jet is 
around this limit, which explains the absence of 
strong shocks. However, the fact that the jet pres- 
sure increases with decreasing distance to the com- 
pact object would generate also a strong shock in 
this simulation if the injection point of the jet in the 
grid were closer to the origin. 

In the case of a jet interacting with a stellar 
wind, the lower limit of the initial jet temperature 
at which a standing shock would appear is reduced 
by an order of magnitude. ^^ This is important, as 
once the jet head is far away, the cocoon will fade 
out and the stellar wind will take its place in the in- 
teraction with the jet, though asymmetrically (only 
from the direction to the primary star), what would 
produce asymmetric standing shocks. 

We have also performed planar jet simulations 
for the weak and mild jet, in which the wind veloc- 
ity is taken into account. Fig. 3 shows the influence 
of a lateral wind on the weak jet. In these maps, 
the asymmetrical standing shocks produced by the 
jet-wind interaction are clearly observed. A second 
simulation for a slab jet with the same physical pa- 
rameters as the mild cylindrical jet shows that a jet 
with power 10^^ erg/s can also be destroyed by such a 
wind.^^ 3-dimensional simulations will be performed 
in order to check this interaction in a more realistic 
setup. 



4. Non-thermal radiation from jet-wind 
interactions 

It is worthy noting the importance of the strong 
shocks predicted by our simulations from the radia- 
tive point of view. On one hand, particle accelera- 
tion is likely to take place in the standing shocks or 
the head of the jet when crossing the denser wind 
regions. The presence of a non thermal population 
of electrons plus the strong photon field produced 
by the primary star will lead to strong IC emission 
(e.g. Ref. 12,13). Wind-jet interactions may be the 
mechanism triggering particle acceleration in sev- 
eral massive X-ray binaries with jets from which 
gamma-rays have been detected. ^^^^ In addition, 
these shocks will possibly increase also the value of 
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the magnetic field (not considered in the simulation 
at this stage) in the postshocked regions to high 
values, allowing for efficient synchrotron radiation 
up to X-rays and beyond (e.g. Ref. 12). Moreover, 
when enough relativistic proton power were avail- 
able, there may be significant gamma-ray produc- 
tion via decay of neutral pions from proton-proton 
interactions due to shock enhancement of the matter 
density, plus other products like secondary pairs and 
even neutrinos (e.g. Ref. 17,18). We note that for a 
10% efficiency for energy transferring from shocks to 
non thermal particles, multiwavelength leptonic lu- 
minosities of ~ 10^^ erg/s, detectable from low (e.g. 
radio) to very-high energies (e.g. TeV), may be easily 
reached in the case of a powerful jet. 
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5. Summary 

We show that, with the reasonable overpressure and 
physical parameters in the jets and stellar winds 
of HMMQ, assuming negligible wind porosities, the 
jets form strong standing shocks due to interac- 
tion with the cocoons formed first, and the stellar 
wind from the primary later, on top of the tran- 
sient reverse/bow shocks. We estimate that stand- 
ing shocks must form within the system region for 
typical jet temperatures (> 10^~^^K). Moreover, we 
have shown that jets with power < 10^^ erg/s could 
be deflected and disrupted by the primary star wind. 
The influence of hydrodynamical instabilities such as 
Kelvin-Helmholtz is to be neglected on the basis of 
our results, as the evolution and dynamics of these 
jets are nonlinear from the start. In this scenario, 
the inertia of the jet compared to that of the wind 
is responsible for the fate of the former. Finally, the 
interaction of the jet head with the wind, and the 
jet sides with the cocoon and the stellar wind, in 
HMMQ, could lead to efficient non-thermal multi- 
wavelength radiation, showing the importance of the 
role of the primary star for high energy emission in 
these sources, which is not just a provider of seed 
photons for IC scattering and photon photon absorp- 
tion, or target nuclei for hadronic processes. 
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Fig. 3. Maps of axial velocity (cm s"-*^), rest mass density 
(g cm~'^) and pressure (dyn) and axial velocity (cm/s) at the 
end of the simulation (tf = 680 s) of a weak jet in slab geome- 
try. The horizontal and vertical coordinates indicate distances 
(cm) to the compact object and the jet axis, respectively. The 
star would be on the vertical axis of the plot, left of the jet, 
at the orbital distance (out of the plot). 



Acknowledgments 
References 

1. Marti, J. M, Miiller, E., Font, J.A., Ibafiez, J.M, and 
Marquina, A. ApJ, 479, 151 (1997). 

2. Saxton, C. J., Bicknell, G. V., Sutherland, R. S., 
Midgley, S. C. J. MNR AS 369, 781 (2005). 

3. Pemcho, M., Marti, J. M., Hanasz, M. A&A 443, 863 
(2005). 

4. Laing, R.A., in Energy Transport in Radio Galaxies 
and Quasars, eds.: P.E. Hardee, A.H. Bridle and J. A. 
Zensus, ASP Conference Series, vol. 100, p. 241 (1996). 

5. Heinz, S. & Sunyaev, R. A&A 390, 751 (2002). 

6. Bordas, P., Paredes. J. M., Bosch- Ramon, V. these 
proceedings (2008). 

7. Zealey, W. J., Dopita, M. A., & Malin, D. F. 1974, 
MNRAS, 192, 731 

8. Corbel, S., Fender, R. P., & Tzioumis, A. K., et al. 
Science, 298, 196 (2002) 

9. Gallo, E., Fender, R., Kaiser, C, et al. 2005, Nature, 
436, 819 

10. Earners, H. J. G. L. M. & Leitherer, C. ApJ 412, 771 
(1993). 

11. Owocki, S. k Cohen, D. ApJ 648, 565 (2006). 

12. Perucho, M. & Bosch- Ramon, V. A&A submitted 



February 2, 2008 9:59 



WSPC - Proceedings Trim Size: llin x 8.5in hepro-draftl 



Log Proper Rest-Mass Density 



1.29e-14 
1.33e-16 
1.38e-18 

3.77e+05 

2.40e+01 

1 .52e-03 



I 



1 10" 
5 10'" 

0^ 
-5 10'° 
-1 10'' 
6.0 10'° 



1 10" 
5 10" 



-5 10 
-1 10" 



9.6 10" 

Log Pressure 




9.6 10" 

Axial velocity 




6.0 10" 



5.1 10' 



9.6 10" 



1.41 10 



1.86 10" 



Fig. 1. Maps of rest mass density (g cm~^), pressure (dyn) and axial velocity (cm/s) at the end of the simulation (tf = 164 s) 
of the mild jet in cylindric geometry. The horizontal and vertical coordinates indicate distances (cm) to the compact object and 
the jet axis, respectively. 
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Fig. 2. Maps of rest mass density (g cm"*^), pressure (dyn) and axial velocity (cm/s) at the end of the simulation (ff = 212 s), 
for a powerful jet. We have included the map of axial velocity in order to illustrate the jump in velocity at the recollimation 
shock. 
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